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In recent years, a number of dynamical density functional theories (DDFTs) have been developed 
for describing the dynamics of the one-body density of both colloidal and atomic fluids. In the 
colloidal case, the particles are assumed to have stochastic equations of motion and theories exist 
for both the case when the particle motion is over-damped and also in the regime where inertial 
qq ■ effects are relevant. In this paper we extend the theory and explore the connections between the 

| microscopic DDFT and the equations of motion from continuum fluid mechanics. In particular, 

starting from the Kramers equation which governs the dynamics of the phase space probability 
£\j ■ distribution function for the system, we show that one may obtain an approximate DDFT that is 

a generalisation of the Euler equation. This DDFT is capable of describing the dynamics of the 
fluid density profile down to the scale of the individual particles. As with previous DDFTs, the 

Q dynamical equations require as input the Helmholtz free energy functional from equilibrium density 
functional theory (DFT). For an equilibrium system, the theory predicts the same fluid one-body 
, density profile as one would obtain from DFT. Making further approximations, we show that the 

theory may be used to obtain the mode coupling theory that is widely used for describing the 
transition from a liquid to a glassy state. 
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I. INTRODUCTION 



o 

To study fluid dynamical phenomena it is often sufficient to consider the fluid as a continuum and ignore the 
fact that it is in reality made up of individual particles. For example, such an approach leads to the Navier-Stokes 
equation, a corner stone in the field of Fluid Dynamics iii2i£ It is only if one is interested in the dynamical behaviour 
on length scales comparable with the size of the individual particles that one must go beyond a continuum theory. 
This is rarely the case for atomic and molecular fluids. However, colloidal suspensions are a particular class of fluids 
where, due to the mesoscopic size of the colloids, this limit is more easily reached. For example, to describe blood 
flow in the capillaries, due to the fact that the diameter of the red blood cells can be similar to that of the capillaries, 
one can not treat the fluid (blood) as a continuum. This regime is also reached in a number of microfluidic devices^ 
and in laser tweezer experiments^ where groups or individual colloids are trapped and moved around by the tweezer. 
Over the last few years a dynamical density functional theory (DDFT) has been developed. This constitutes a 
t—\ , microscopic theory for the fluid dynamics of such colloidal fluids. The starting point was work by Marconi and 
7-H ' Tarazona^i in which they assumed that the colloids can be modelled as Brownian particles with stochastic equations 
of motion, thus neglecting hydrodynamic interactions between the colloids. Newton's equations of motion for a system 
■ of N identical Brownian particles of mass m, specified by the set of particle position coordinates r N = {ri , r2, • • • , rjy} 
00 ' and momenta p N = {pi, P2, ■ • ■ , Pat}, are£ 

o 



dvj _ Pi 

dt m ' 

^=- 7Pi + X i (r i )+G i (t), (1) 

where the forces due to the solvent on the colloidal particles are modelled solely by a viscous drag term, — 7Pi, where 7 
is a friction coefficient, and by Gj(i) = (£f (i), £f (t), £f(t)), a stochastic white noise term with the property = 

and (t')) = 2"fmkBT5ij5°"' 'S(t — t'), where fc# is Boltzmann's constant and T is the temperature. X, is the 

sum of the force on particle i due to any external potentials and the forces due to interactions with the other particles 
in the system, and is given by Xj = — V ri V(r N , t), where the potential energy V(r N , t) is assumed to be of the form 

V(r N , t) = V^{v u t) + i £. j v 2 (r, t (t), Tj (t)) 

+gE i , i , t "3(ri(i),r i (f),r fc (i))+- , (2) 

where V ext (r,t) is the external one-body potential, V2 is the pair interaction between the particles, V3 is the triplet 
interaction (which is often assumed to be zero) and there may also be other higher body terms which are denoted 
by • • • in Eq. ([2]). We should emphasise that the hydrodynamic interactions between the colloidal particles, which can 
be vital for a quantitative description of the dynamical processes in dense suspensions, have been neglected in Eq. 
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([T]). To include these effects, the term in Eq. (Q} describing the force exerted by the solvent on colloid i, F^ 1 = — 7Pi, 
should be replaced by F^ 1 = — Y^,jLi ^ij( rN ) ' Pj> where Yjj(r N ) is a microscopic friction matrix that depends on the 
coordinates of all the particles.— 

There are two limits one may consider in Eq. ([1]). The first is the limit when 7 is small. In the limit 7 — > 0, both the 
viscous drag term and the stochastic noise term in Eqs. |T]) disappear (recall that they are related to one another by 
a fluctuation-dissipation relation) and Eq. {1} reduces to a (deterministic) set of Newton's equations for the motion 
for a system of particles that are not immersed in a solvent: 

m^=X,(r«). (3) 

These are the equations of motion for an atomic or molecular fluid. 

The other limit to consider is when the friction coefficient 7 is large (overdamped dynamics). In this case the term 
dpi/dt in Eq. ([T]) is small and one may neglect it, so that Eq. ([T]) reduces to the following set of equations of motion 
for the colloids: 

dr- 

7 m^=X i (r < ) + G i (i). (4) 

For a system with these underlying equations of motion, Marconi and Tarazona obtained the following equation of 
motion for the one-body density p(r, t) of the colloidal particles&iiii 



dp(r,t) _ 1 y 
dt 7m 



5p(r,t) 



(5) 



by making the assumption that the two-body correlations in the non-equilibrium fluid can be approximated by those 
of an equilibrium fluid with the same one-body density profiled In Eq. ((SJ, F[p] is the equilibrium Helmholtz free 
energy functional from the well-established classical density functional theory . 13 ' 14 It is the fact that the theory builds 
upon the equilibrium free energy functional that makes the theory so appealing. It means that so long as one has a 
reliable expression for F[p] (and many such approximations are known - see, for example Refs. [l3l [Til and references 
therein) then Eq. ([5]) at least guarantees the correct equilibrium profile^^ and in practise, it often proves to be 
reliable out of equilibrium as well^^^ii&ilii&i&S&Z 1 

The main reason that one can obtain such a relatively simple expression as Eq. ([5]) for the equation governing the 
fluid dynamics is because the colloids are suspended in a solvent, which acts as a heat bath, so there are no thermal 
gradients and one does not explicitly have to take into account conservation of energy. Neither does one have to 
account for momentum currents, due to the overdamped dynamics. Recall that such considerations are required to 
obtain reliable dynamical equations for atomic or molecular fluids such as the Navier-Stokes equation^ However, even 
for colloidal fluids one should take account of momentum currents, since strictly speaking, Eq. {1} provides a better 
account than Eq. Q of the particle dynamics. It was to address this issue that Marconi and Tarazona published a 
further paper, 22 in which they derived, using a multiple time scale analysis, a DDFT that incorporates inertial effects. 
In the limit where 7 is large, their theory reduces to their original theory, 6 i.e. where the dynamics is governed by 
Eq. . Their analysis generates a hierarchy of coupled equations which in principle must be solved self consistently. 
They also argue that so long as 7 is not too small, then one can describe the fluid dynamics using just the truncation 
of the expansion at second orderj 22 i 23 

In this paper we take a different route and derive an alternative DDFT, applicable to the same kinds of systems. 
We do not use the multiple time scale analysis that Marconi and Tarazona employed^ Instead, starting from the 
Kramers equation for the iV-particle phase space probability density distribution function of the system f^ N \ we 
derive an equation (Eq. (|16[1 ). which describes the dynamics of the momentum currents j in the fluid of Brownian 
particles. In this exact equation, we make two approximations: (i) We assume that one can approximate the two 
body spatial correlations in the non-equilibrium fluid by those of an equilibrium fluid with the same one body density 
profile - this is also the approximation used to obtain Eq. ([5]) . (ii) We assume a local equilibrium Maxwell-Boltzmann 
form for the one particle phase space probability density distribution function f^'. From these two approximations 
we are able to derive a DDFT that takes the form of a generalised Euler equation. The Euler equation may be 
obtained from the Navier-Stokes equation by setting the shear and the bulk viscosity to zero. The same is true for 
the present system: assumption (ii) above is equivalent to setting these viscosities to zero, and following Kreuzer— we 
see that by going beyond the local equilibrium approximation for J*- 1 -*, we may obtain a DDFT that takes the form 
of a generalised Navier-Stokes equation. These two equations from fluid mechanics are very well known, and provide 
the starting point for the description of many fluid dynamical phenomena. A large body of knowledge has been built 
up in the literature concerning solutions of the Navier-Stokes and other such equations from fluid dynamics The 
contribution of the present paper is to show how to connect these continuum theories with the fully microscopic 
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DDFT, and to indicate how one may incorporate information about the microscopic structure and correlations of the 
fluid into the continuum theories. 

To solve both the generalised Euler and generalised Navier-Stokes DDFTs mentioned above, requires explicitly 
keeping track of both the fluid density p(r, t) and the local fluid velocity v(r, t). However, by neglecting certain terms 
in the dynamical equations, we show that one may obtain a DDFT that only explicitly depends of the fluid density 
p(r,t). From this DDFT, following the approach given in Ref. H3, we are able to obtain a mode coupling theory 
(MCT) for the density fluctuation correlation function, that is of the standard MCT form. 

This paper is laid out as follows: In Sec. HH we derive a DDFT for the average one-body density p(r, t) of a fluid of 
colloidal particles, whose equations of motion are given by Eq. fl}. In Sec. IIIII we use our DDFT to obtain the MCT. 
Finally, in Sec. II VI we summarise and draw some conclusions. 



II. EQUATIONS OF MOTION 

For a system of N identical particles, whose equations of motion are given by Eq. (JTJ) , the time evolution of the phase 
space probability density function, f( N '(r N ,p N ,t), which gives the probability of the system being in a particular 
configuration (r , p-^) at time t, is governed by the Kramers (Fokker-Planck) equation^ 



Qf(N) 



N N N N 



()l - ^EprV n /W+^X r V Pi /W= 7 ^V Pj . Pi /W + 7*T^V^W. (6) 

i—1 i—1 i—1 i—1 

Note that in the limit 7 — > this reduces to the Liouville equation^ which is the equation governing the time 
evolution of f( N ' when the system equations of motion are given by Eq. (J3J) . We may define a set of reduced phase 
space distribution functions^ 



/W(r«,PV) = (^yy / dr<"-*> J dp<"-»>/W(r", p", t). (7) 

Integrating over the Kramers equation ([6]), we obtain the following dynamical equation for the one particle reduced 
phase space distribution function /W; 

§i + m ' Vri +Fft ' Vpi ") ^ (1) ( ri ' Pl '*) = i [ y pi ' Pl +mk B TV 2 pi ] /«(n, Pl ,i) 



J dr 2 J dp 2 F 12 - V Pl / (2) (r 1 ,p 1 ,r 2 ,p 2 ,t) + -- - , (8) 



where is the two-particle distribution function, Y\ xt — — V ri V ext (ri), F i2 = — V ri u 2 (ri — r 2 ) and • ■ • contains 
contributions from three-body and higher-body interactions. We have also assumed that all boundary terms are zero 
- i.e. that /(™)(r n , p n , t) is zero when any of the components of and pi — > ±00. Note that when 7 = 0, Eq. 
((SJ) is simply the first equation in the Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy^ The following 
manipulations are similar to those in Refs.Q] andl24l although the starting point, Eq. 1(5)), is different. In Eq. (jHJ) we 
may integrate over the momentum pi, in order to obtain the continuity equation: 



^ + V ri .j = 0, (9) 



where the one-body (number) density 



and 



p(ri,t) = J dpi/W(n >P i,i), (10) 

j(r 1) t)= f d Pl 5i/«(ri lPlj i) (11) 
J m 

is the current. Taking Eq. |8[). multiplying through by pi/m and then integrating over the resultant equation with 
respect to P i, yields the following: 
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ff(ri,f) 
dt 



J m z m 



-— { dr 2 Fi 2 p (a) (ri > r2,t)-— / dr 2 / dr 3 F 123 p (3) (ri, r 2 , r, 3 , *) + • • • = 0, (12) 
m J m J J 

where F123 = — V ri U3(ri, r 2 , r 3 ) is the three-body force on particle 1 due to particles 2 and 3, 

p^(r ll r 2 ,t)= J dpi /" dp 2 /( 2 )(r 1 ,pi,r 2 ,p 2 ,t) (13) 

is the two body density distribution function and 

p^(n,r 2 ,r 3 ,t) = J dpi J dp 2 J dp 3 / (3) (ri,pi,r 2 ,p 2 ,p3,t) (14) 

is the three body density distribution function. Note that pi (8 Pi in Eq. (|12p denotes a tensor product (dyadic). 
Eq. ([9]), the continuity equation, is simply a statement of conservation of mass and Eq. (|12p is a momentum balance 
equation. These two equations are just Eqs. (19) and (20) in Ref. HH, and it is from here that we take a different 
approach to the authors of Ref. in developing the theory. 

At equilibrium, the one particle distribution function takes the Maxwell-Boltzmann form: 14 

f"H r .P>= * r > «v(-^). (15) 



(2nmk B T) 3 / 2 * \ 2mk B T 

Thus, at equilibrium, we find that the integral J dpi(pi ® Pi)/^ 1 - 1 = Tnk B T p{rx)\, where 1 denotes the 3x3 unit 
matrix. Noting this result, we may recast Eq. (fT2|) as follows 

+ 7 j(n, t) + A(n,t) + ^Vp(r 1( t) + -p(r u t)V^ 4 (ri, t) 
at mm 

+ — / dr 2/ 9 (2) (ri,r 2 ,i)Vv 2 (ri - r 2 ) 
m 7 

/" dr 2 y dr 3/9 (3) (ri,r 2 ,r3,t)Vv 3 (ri,r 2 ,r 3 ) + --- = 0, (16) 

where 

A(r,t)=V./dpte-*^l)/«(r, P) t). (17) 
J \ m m J 

From equipartition, it can easily be seen that at equilibrium A(r, t) = 0. Note that in the limit 7 — ► 0, i.e. when the 
solvent is absent, Eq. ()16|) is still correct, even though the underlying equations of motion |T]) become deterministic. 

So far, we have made no approximations and together Eqs. (fit?)) and (fT7|) are exact. It is at this stage in the 
derivation that we make two approximations. The first is to assume that we may approximate the two-body spatial 
correlations in the non-equilibrium fluid by those of an equilibrium fluid with the same one body density profile. 
This is done by assuming that one can apply the following sum rule (which is exact for the equilibrium fluid) as an 
approximation for the sixth and seventh terms on the left hand side of Eq. (|16p >^'^ 

-k B Tp(r 1 )Vc^(r 1 )= J dr 2j0 ( 2 )(n, r 2 )V ri u 2 (ri, r 2 ) 

+ y dr 2 y dr 3 p (3) (n,r 2 ,r 3 )V ri w 3 (ri,r 2 ,r 3 ) + -- - , (18) 

where c^(r) is the one body direct correlation function and is equal to the functional derivative of the excess part of 
the Helmholtz free energy functional i 13 ' 14 

c(1)(r) _ (19) 

where = l/k B T. 
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Assuming that we may apply Eq. (| 18|) for the non-equilibrium fluid should be reliable for the case when the 
particles interact via potentials that are slowly varyingj2i22 However, in cases such as a system of hard spheres, 
where the collision dynamics are somewhat different than in systems of particles interacting via potentials that vary 
continuously, this approximation is less reliable, and one should approximate this quantity using a term involving a 
binary collision operator, along the lines described in Refs. |22| andl23l 

Making this approximation in Eq. (|16j) . and also noting that Vp = pVln(p), we obtain: 

+ 7 j(r, t) + A(r, t) + Ip (r> t)V 5 -^^ = 0, (20) 
ot m op(r,t) 

where the Helmholtz free energy functiona l 13 ! 14 

F[p] = k B T J drp(r)[ln A 3 p(r) - 1] + F ex [p] + J dvV ext (r)p(v). (21) 

The first term on the right hand side is the ideal-gas contribution to the free energy, A is the thermal de Brogle 
wavelength and -F e a:[p] is the excess (over ideal-gas) contribution due to interactions between the particles. 

We now make a second approximation and we assume that we can make a 'local-equilibrium' Maxwell-Boltzmann 
approximation for the one-particle distribution function 1 -^ [c.f. Eq. (|15[)]: 

(2nmkB-i ) 1 V ImRB-L J 

where p(r, t) — mv(r,t), and v(r, t) is the average local velocity of the particles. If one substitutes Eq. (f2"2")) into 
Eq. (fTTT) then one obtains the following expression for the current: 

ji. e .(r,t) = p(r,t)p(r,t)/m 

= p(r,t)v(r,t). (23) 

Using this approximation for the current, the continuity equation (|9]) becomes: 

g + V-(pv)=0. (24) 



On substituting Eq. ([22]) into Eq. ([T7|. we obtain the result 

>M) 



A u Xr,t) = V 



p(r,*)®p(r,t) 



(25) 



This, combined with Eq. (|2H)) . gives: 



9(PV) + 7 pv + V • [pv ® v) = -ipV^l. (26) 



9i ' m 5p 

This expression may be simplified by noting the following results 



V • (pv ® v) = p[V • (v ® v)] + (Vp) • (v ® v) 

= p[v(V • v) + v.Vv] + v(v • Vp) (27) 



where Vv is the tensor derivative, and 



c>(pv) <9v 9p 

= P ~dt +V dt 
d~v 

= p— - pv(V ■ v) - v(v ■ Vp) (28) 

where we have used Eq. (|2"4")l to obtain the second line in Eq. (|28p. We now use Eqs. (|2T|) and (|28p to simplify Eq. 
(US}, giving: 

p( ^+v-Vvj + 7 pv = --pV^p. (29) 
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Dividing through by the fluid density we obtain: 



Dv 1 SF\p] , . 

— + 7 v = V — ^. 30 

Dt m op 



where 

D d 



Dt-dfW (31) 



is the substantive derivative. Taken together, Eqs. (|24)) and (|29)) are one of the main results of this section of the 
paper. Given a suitable approximation for the free energy functional -Fe^fp] in (f2Tjl . Eqs. (|24|) and (|29p may then be 
simultaneously solved (with suitable boundary conditions) for the dynamics of the system. 

It is worth noting at this point that when 7 is large, the first term on the left hand side of Eq. (|30|) is negligible 
with respect to the second term on the left hand side, giving 

1 SFlp] , x 

7v^ V— ^. 32 

m op 

On substituting this expression into Eq. (f2"4"| , we obtain Eq. ([5]), the original DDFT of Marconi and Tarazona, 6,7 
In the opposite limit, when 7^0, Eq. ([25]) reduces to a generalised Euler equation: 

„,„(£ +v . V v)=-,vfd (33) 

(recall that here we use the notation that p(r,t) is the colloidal fluid number density). If we make a local density 
approximation (LDA) for the Helmholtz free energy functional (pTTl) : 



F[p] = J dxf(p(r,t)) + J drV ext (r)p(r,t), (34) 
where f(p) = kgTp^lnA 3 p — 1) + f ex (p) is the intrinsic Helmholtz free energy density, then from Eq. (|3"3"|) we obtain 



mp[ — +v-Vvj = -pV^-pV^* (35) 



which may be rewritten in the following way: 

mp(^+v-Vv) = -Vp- P X7V ext , (36) 

where p = — f + pp is the local pressure and p, = df /dp. At equilibrium, p is the chemical potential. This equation 
is the Euler equation from Fluid Dynamics!^ 

We are now in a position to see the consequence of making the local equilibrium approximation (|22p for the one 
particle distribution function / W (r, p, t) in Eq. ([20]) . This approximation leads to effectively setting the shear viscosity 
7] = (recall that one obtains the Euler equation by setting r\ = in the Navier-Stokes equation) . To go beyond this 
approximation, one may follow Kreuzer— and assume that the distribution function /^'(r, p, t) can be expanded as a 
Taylor series, as follows: 



/« (r, p, t) = /^(r, p, t) + ai (\p - mv|) 



(p — mv) ■ A • (p — 772v ) 1 

(p — 777V) 2 3 



(37) 



where the function a\ may also be a function of position and A is the symmetric rate of strain tensor, whose components 
are 

With this approximation in Eq. (|20p and further assuming that the one-body density of the colloids is a constant 
(i.e. that it is an incompressible fluid), we obtain an additional term in Eq. l|2T)|) that is ~ — 77^^V 2 v, where rj^ K > is 
the kinetic energy contribution to 77, the coefficient of shear viscosity, and is formally given by an integral over the 
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distribution function ai(p). 1 Similar considerations at the two body level in the interaction terms in Eq. (|16p give an 
addition contribution ~ t^t/^V 2 v i The resulting dynamical equation 

mp + v ■ Vv) + 7mpv = -pV + r?V 2 v, (39) 



at 1 J ' ' ' ' s P 

where rj = + rj^ v \ is a generalisation of the Navier-Stokes equation for an incompressible fluid. For colloidal 
fluids, we believe that taking Eq. (|3"9"|) together with the continuity equation (f2"4")) may provide a basis for reliably 
describing the fluid dynamics, even when the density is not a constant (i.e. when the fluid is compressible), ft should 
be noted, however, that although one may obtain formal expressions for the quantities rf K ^ and r/ v ) (see Ref. 0] for 
further details), in practice it may be necessary to obtain this quantity by other means. 

Finally in this section, we remind the reader that at equilibrium, where v = and dp/ dt = all of the dynamical 
equations derived above yield the following expression: 

V^M = 0, (40) 
dp 



which may be integrated to give 

SF[p] 
Sp 



M, (41) 



where p. is the chemical potential. This equation for the one-body density profile p(r, t) is exact i 13 ! 14 Thus, given a 
reliable approximation for F[p], all of the theories presented above exhibit the key feature that at equilibrium they 
yield the correct fluid one-body density profile. 



III. DYNAMICS OF THE ONE-BODY DENSITY AND CONNECTIONS TO MCT 



One of the most appealing features of the original DDFT of Marconi and Tarazona, 6 Eq. ([5|), is that it gives a 
description of the fluid dynamics solely in terms of the one body density p, as opposed to the description given in the 
previous section involving both the density p and the fluid velocity v. We now show how to obtain an approximate 
theory that just involves the fluid density p. We take the divergence of Eq. (fll))) and use Eq. © to eliminate the 
terms involving dj/dt and j to obtain: 



d 2 p(ri,t) 

dt 2 



dt 



— V 2 p(r, t) + - V ■ [p{r, t)X7V ext (r, t) 
m 



+— V- / dr'p (2) (r,r',t)Vu 2 (r-r') 
m J 

Jdr' J drV 3) (ri,r 2 ,r,3,i)V« 3 (r,r',r") + -- - . 



(42) 



So far, no approximations have been made and, taken together, Eqs. (|4"2"|) and (fT7|) are exact. If we make the 
approximation used earlier, where we assume that we may use Eq. (|18p as an approximation for the terms involving 
the non-equilibrium fluid pair, triplet and higher correlation functions, then we obtain: 



d 2 p dp 
dt 2+7 di 



V • A(r,t) = — V 



p(r,t)V 



SF[p(v,t)} 
Sp(r,t) 



(43) 



Since A(r, t) ~ j <8> j, when p/V2mfcsT is small, A(r,i) will also be small and we expect that in such circumstances 
this term may be either neglected or its influence incorporated into a renormalised 7* - i.e. 7 is replaced by the 
renormalised 7* = ^, where v = fcsT/mD and D is the self diffusion coefficient^ Neglecting the term involving 
A(r,i) in Eq. (JUJ) gives 



dt 2 dt m 



p(r,t)V 



SF[p(r,t)} 
8p(r,t) 



(44) 



The appealing feature of this equation is that it gives a theory for the fluid dynamics solely in terms of the fluid 
density p. 
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One circumstance where the term involving A(r, t) may be neglected is when the fluid density is high and the 
system is not too far from equilibrium (i.e. when the average local velocity v is small). This situation was explored 
in Ref. for the case when 7 = 0, where it was shown that one may derive the MCT (a theory for a density 
fluctuation correlation function), starting from Eq. (1441) . The argument presented in Ref. |24| is entirely applicable 
to the present case. We do not repeat the full argument here, merely reminding the reader of the salient points. 
One key issue is that we now interpret the one body density p as a coarse grained density field. For example, we 
may follow Refs. [H andl24l and define the temporally coarse grained density p(r, t) — dtK(t — t')p(r, t'), where 

p(r, t) — J2iLi <K r — r i W) i s the density operator (recall that r^t) is the location of the ith particle at time t), and 
K(t) is a normalised function of finite support. The precise shape of K(t) defines the degree of coarse graining. 
The coarse grained density p(r, t) exhibits thermal fluctuations and the amplitude of these depend on the extent of 
the coarse graining*! 2 - The dynamical equations for the coarse grained density depend on a coarse grained two-body 
density distribution function p^(r, r', t) = dtK(t — t')p(r, t')p(r' , t') (and in the case when V3 ^ it also depends 
on a similarly defined three-body distribution function). If we then assume that we have coarse grained sufficiently 
that Eq. (p~5|) still holds true for this coarse grained p^ (r, r', i), then following the argument presented above, we 
obtain Eq. (|44[) as the equation governing the time evolution of the coarse grained one-body density profile. 

The exact excess Helmholtz free energy functional (HU, that is required as input to Eq. (j4~4"|) , is unknown. However, 
we may approximate this quantity by Taylor expanding in powers of Sp(r, t) — p(r, t) — p/,, where pb is the bulk density 
of the uniform fluid. Truncating the expansion beyond terms of 0(Sp 2 ) we obtain; 13 i 24 

F ex [p(r,t)} =F ex [ Pb ] -c«(oc) J dr5p(r,t) - ^ / dr J dr'Sp(r,t)Sp(r' ,t)c^(r - r'), (45) 

where 

(2)/ i\ a S 2 F ex [p] 

C ^ r) = P 5p(r)5p(r>) (46) 

is the pair direct correlation function i 13 i 14 On substituting Eq. (|4"5|) into Eq. ([4^]) and then Fourier transforming, we 
obtain)^ 

k 2 Ohk 2 IIP 

/5k (t) +7Pk(i) = —s-pk(t) + ^— c kPk (t) + -^7^3 / dk'k • kW(*)ck'/Sk-k'(*), (47) 
pm pm pra (ztt)* J 

where Ck is the Fourier transform of c' 2 '(r) and pk(t) is the Fourier transform of Sp(r,t). We multiply through in 
Eq. (|4T)) by p_k(0), and then average over the ensemble of initial configurations of the density field. Making the 
approximation proposed in Ref. [24] (see also Ref. I25T) for the quantity (i?k(^)-R-k(0)), where 

Mt) = } f dk'k • kW(t)ck'flk-k'(*), (48) 

we obtain the MCT equation 
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M*) + 70k(i) + VIM*) = - / dt'mk(t')0k(t - 0. ( 49 ) 







for the normalised density fluctuation correlation function 

(Pk(O)p-k(O)) 

where il 2 , — k 2 / fimS-k and Su = [1 — pbC^]^ 1 is the static structure factor. Making the usual MCT approximation of 
factorising four-point correlation functions into products of two-point correlation functions one obtains the standard 
MCT expression for the memory function! 14 i 24 ' 26 

™k(i) = o^fa / dk ' I k ' k ' c k' + k • (k- k')c k - k '] 2 Sk/S k _ k ^(i)^k-k'(t). (51) 
2(271")^ k z ■m J 

The MCT obtained above for under-damped colloidal particles formally has the same structure as the standard MCT 
that is used to elucidate the properties of the glass transition in atomic/molecular fluid a 14 i 26 - i.e. systems with 
deterministic Newtonian dynamics - see Eq. ([3]). What the above analysis shows is that results obtained for atomic 
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and molecular ffuid a 1 i 26 are also relevant to under-damped colloidal suspension. Since the dynamics in dense glassy 
systems is strongly influenced by how collective density fluctuations decay, it is perhaps not surprising that the decay of 
(f>k(t) should be governed by equations having a structure which do not strongly depend on the underlying equations 
of motion for the individual particles. Note also that for overdamped colloids (large 7), we may neglect the term 
d 2 p/dt 2 from Eq. (|44"]) . to obtain Eq. ([5]). On following the argument presented above, starting from Eq. ([5]), we 
obtain a MCT equation the same as Eq. (|4"9")) , but with the term (j>k(t) omitted from the left hand side. This MCT 
for overdamped colloids was originally obtained by Szamel and Lowen^ 



IV. SUMMARY AND CONCLUDING REMARKS 



In this paper we have derived a number of approximate DDFTs, Eqs. ([29]). (|39p and (|44[) , for systems of Brownian 
particles where inertial effects are relevant. These theories are also relevant to atomic and molecular fluids, which 
correspond to the limit 7 — > in these equations. To derive these DDFTs we made two approximations. The 
approximation that pertains to all three cases is to assume that one may use Eq. (|18[) . which strictly only applies at 
equilibrium. In using this, we effectively assume that the two-body spatial correlations in the non-equilibrium fluid 
are equal to those in an equilibrium fluid with the same one body density profile. As discussed in Refs. and H^, this 
is essentially a mean field approximation that is reliable only in the high friction limit or for systems interacting via 
pair potentials that are continuous and differentiable. In Refs. [22| and [23 the authors consider the particular case of a 
system of hard spheres. The collision dynamics of such a system is somewhat different to that in systems of particles 
interacting via potentials that vary continuously. They show that to describe the hard-sphere collision effects, one 
must approximate the term on the right hand side of Eq. (| 1 2jl involving the two-body distribution function by a term 
that contains a binary collision operator— (the three-body potential is zero) . Such an approach may also prove to be 
useful for extending the theory presented here. 

The second approximation that was made here was to assume a particular form for the phase space probability 
distribution function - see Eqs. (|2"2"|) and ([57)) . Making approximations at this level are justified for systems that 
are not too far from equilibrium. However, for very strongly driven systems, one must expect to have to go beyond 
such approximations to obtain a reliable description of the fluid dynamics. 

The approach presented here provides a fully microscopic basis for the equations used to consider systems with 
diffuse interfaces in fluid mechanics. Either one-component (gas-liquid) systems may be considered,— or the present 
results may easily be extended to consider a two component fluid, in order to study the dynamics of a system 
exhibiting fluid-fluid phase separationj 30 i 31 i 32 These approaches may be obtained by assuming a gradient expansion 
of the Helmholtz free energ y 13 ' 1 in the present theory. As mentioned in the introduction, and we emphasise again 
here: the DDFTs presented here are microscopic generalisations of the Euler and Navier-Stokes equations. Much is 
known in the literature about the solutions of these equations^i 3 - One such method worth mentioning, due to having 
some underlying connections to the present approach, is the Lattice-Boltzmann methodj 37 ' 38 The present work shows 
how to build upon this knowledge and to incorporate into these theories information about the microscopic structure 
and correlations in the fluid via the Helmholtz free energy functional. 

A further extension of the theory presented in this paper is to use the DDFT to obtain a MCT theory that is 
applicable for studying the glass transition in colloidal suspensions in which inertial effects in the particle dynamics 
are important. To derive the MCT, we assume that we may neglect the term V • A(r, t) in Eq. (J33]). This term is 
significant when currents are large (i.e. when the condition v <C ^Jk^Tjm, no longer holds) and so makes the theory 
unreliable in such cases. However, for dense fluids that are near to equilibrium or with a low average velocity v, one 
should find that the DDFT in Eq. (|44[) is relevant. We should also mention that stochastic dynamical equations, 
that are of a similar structure to the dynamical equations that we have obtained here (Eqs. ([29]) . ([39]) ). have been 
developed to study the dynamics in glassy systems . 33 i 34 i 35 i 36 Within the present DDFT framework, these stochastic 
analogues of Eqs. (f29|) and (l39|) may be viewed as having been obtained by considering the dynamics of a coarse 
grained density field, along the lines presented in Ref. 

We should also make a few further comments concerning Eq. (|44[) . Equations of this general form have appeared 
before in the literature in a number of different contexts. When the free energy functional F[p] in Eq. ([44]) is set to 
be simply that of an ideal gas (i.e. when F ea; [p] = 0) and we set the external potential V ext (r) = in Eq. ([2"Tj) . then 
we obtain 

d 2 p dp k B T 2 

W +1 m=— Vp - (52) 

This equation is telegrapher's equation and has numerous applications. Equations of this general form have been 
used for example in the description of heat waves^ electric currents 40 and nuclear collision dynamics- 41 When F[p] 
is replaced by the Ginzburg-Landau free energy functional, a generalised Cahn-Hilliard equation is obtained. Such 
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an equation was proposed by Galenko and co-workers as a phase-field model for binary alloy a ' i ' and also by 
Koide et. a?j 46 i 47 as a means of incorporating memory: if one assumes overdamped particle dynamics (Eq. and 
a noise field Gj(i) with memory - i.e. where (£f ^ %<5 a!/ exp(— \t — t'\/r), then following the argument of 

Koide et. al, one obtains an equation of the same form as Eq. (j4"4")) . 

To conclude, we recall that the present theory is for systems with underlying equations of motion given by Eq. |T]). 
Therefore, the DDFT will not be reliable for systems where the equations of motion are not well- modelled by Eq. |lj, 
such as in colloidal fluids where hydrodynamic interactions between the colloids are significant. One may include the 
hydrodynamic interactions by extending the original Marconi- Tarazona DDFT ^ to include additional terms that 
describe the hydrodynamic interactions at the Rotne-Prager level4& Alternatively, one may perhaps be able to use 
the present DDFTs to incorporate hydrodynamic effects by treating such systems as two component mixtures. The 
hydrodynamic interactions between the colloids would enter the treatment via the density field of the second (solvent) 
species. 

I gratefully acknowledge helpful discussions with Lubor Frastia, Markus Rauscher and Uwc Thiclc and thank RCUK 
for financial support. 
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